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We review the Rein-Sehgal model and criticize its use for low energy neutrino induced coherent 
pion production. We have studied the validity of the main approximations implicit in that model, 

' trying to compare with physical observables when that is possible and with microscopical calcula- 

tions. Next, we have tried to elaborate a new improved model by removing the more problematic 

■ approximations, while keeping the model still reasonably simple. Last, we have discussed the limi- 

tations intrinsic to any approach based on the partial conservation of the axial current hypothesis. 
In particular, we have shown the inability of such models to determine the angular distribution of 
the outgoing pion with respect to the direction of the incoming neutrino, except for the (f — Q 
kinematical point. 

! PACS numbers: 25.30.Pt,13.15.+g,12.15.-y 

! I. INTRODUCTION 

The theoretical modelling of neutrino-induced coherent pion production in nuclei in the low (~ 1 GeV) energy 
region is a quite compHcated task. The reasons are twofold. First ^ur limited knowledge of the weak pion production 
on nucleons at these energies, (see for instance Refs. P, i, i, H E i, 0, i, H), that relies upon conflicting and low 
statistics bubble chamber experimental data (lol . [llj . Second, the need of a quantum mechanical treatment of the 
multiple scattering involving the strong interaction of pions and nucleons that can hardly be accommodated in typical 
MonteCarlo approaches. Many different models have been proposed. Some of them are based on the assumption that 
coherent pion production is dominated by the divergence of the axial current and the use of the Partial Conservation 
of the Axial Current (PCAC) hypothesis [11, [11, [Ti, [11, [11, [11 [11, . Other works develop a microscopical model 
for both pion production and distortion assuming that the pion nuclear interaction is dominated by the A(1232) 
resonance [l^, [2l|, [H, [H, [13, [1^ • Thus, they can only be applied at relatively low neutrino energies. Despite the large 
effort dedicated to this process the situation is hardly satisfactory and large discrepancies are found among different 
theoretical predictions, both for total and differential cross sections. 

Besides, the understanding of these processes is very important for the analysis of neutrino oscillation experiments. 
For instance charged current (CC) coherent pion production is one of the candidate processes responsible for the deficit 
found by the K2K collaboration in the forward scattering events, which seriously limited the prediction accuracy of the 
neutrino energy spectrum at the far detector |26l] . Similarly, neutral current (NC) 7r° production is one of the largest 
background sources in muon neutrino experiments, in particular for electron neutrino appearance measurements 
' [27l[28|. because it can mimic electron events when one of the 7r° decay photons is not detected. A significant part of 
the tt'^'s in the forward direction could come from coherent production. 

Due to its success in the description of CC coherent pion production data at high energies and also its simplicity and 
easy implementation, the Rein-Sehgal (RS) model of Ref. fl^ has been widely used by experimental collaborations 
in their analyses j26. , .28. , .29i] . even at quite low neutrino energies clearly beyond the scope of the original paper. We 
should remark that some of the RS model basic approximations are better suited for high energies and heavy nuclei, 
as this implies more forward peaked cross sections. In fact, the minimum neutrino energy discussed in Ref. [12| was 2 
GeV whereas the model has been used for the analysis of the coherent process at neutrino energies down to the pion 
production threshold. 

When compared to other approaches, which are more reliable for low energy neutrinos, the differences clearly show 
up in the total production rate (see e.g. Fig. 2 of Ref. [2§|) and also in the pion angular or energy distributions (see 
Fig. 8 of Ref. ^5]^). Some of the discrepancies, like the wider angular distributions of the pions predicted by the 



^ The curves labeled as MiniBooNE Coll. correspond to the MiniBooNE implementation of the RS model 
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RS model, can be easily traced back to the approximations used in its derivation and they could lead to important 
consequences in the determination of the size of the pion yield and/or the ratio of coherent to non-coherent pion 
production for low energy incident neutrinos. 

Given the deficiencies of the RS model for low neutrino energies, it seems clear the need of the use of alternative 
approaches. However, up to our knowledge, the currently available models are of limited value for that purpose. On 
the one hand, the PCAC based calculations share many of the limitations of the RS model that will be discussed 
below. On the other hand, the microscopical models based on the A dominance are expected to be reliable only for 
low energy neutrinos 1 GeV). 

Some important steps have been given in the direction of developing a model for pion production on the nucleon 
beyond the A resonance, e.g. Ref. [30|, but we are still in need of more experimental data (cross sections) that could 
better constrain the little known transition form factors. This is a long term project that will require a collaborative 
effort of both experimental and theoretical groups. 

In the meantime, and lacking a microscopical approach for the whole energy region of interest, it looks worthwhile 
to try to modify the RS model to extend its applicability towards lower neutrino energies. Some recent works in this 
line include Refs. [H, [3l| which already correct some of the known RS problems. In this work, we will further explore 
this possibility. First, we will carefully study the validity of the main approximations implicit in the RS model, trying 
to compare with physical observables when that is possible and with microscopical calculations. Second, we will try 
to elaborate a new improved model by removing the more problematic approximations, while keeping the model still 
reasonably simple. Last, we will consider the limitations intrinsic to any PCAC based approach. 

The paper is organized as follows. After some basic formalism in Sect. [Tll we review the RS model in Sect. IIIII 
In Sect. IIVI we present our improved model. In Sect. |V] we confront an extension of both models with low energy 
pion- nucleus scattering data. In Sect. lVli we show the results of both models for low energy neutrino- induced coherent 
pion production and we also compare them with a microscopical calculation that we consider reliable in that energy 
region. Finally, we summarize the paper in Sect. IVIII 



II. NC NEUTRINO COHERENT PION PRODUCTION STRUCTURE FUNCTIONS 

We will focus on NC tt*^ coherent production off a nucleus in its ground state Afgs , 

iyiik)+Mgs^Mk')+Mgs+7r°ik^). (1) 

The modifications required for the CC case are straightforward. The process starts with a weak pion (tt^) production 
followed by the strong distortion of the pion in its way out of the nucleus. The nucleus is left in its ground state 
unlike the case for incoherent production where the nucleus is either broken or left in some excited state. 

Defining the four momentum transfer q ~ k—k' and taking q and kxk 'along the positive z— and y— axis respectively, 
one can write the differential cross section with respect to the outgoing neutrino variables, the Mandelstam variable 
t = {q — fc^)^ and the pion azimuthal angle, (f>k„q, in the LAB system as: 

dE'dn{k')dtd(f>k„q \k\ 167r2 ' 
with E' the energy of the final neutrino and G the Fermi constant. The leptonic tensor is given by: 

Lfia = k^kfr + k^k^ H ^5^"" ^ ^^fiaapk k^ , (3) 

and it is orthogonal to the transferred four momentum q^. In our convention, we take £0123 = +1 and the metric 
g^iv _ (^_j^^_^_^_y 'pj^g hadronic part, 7i, is determined by the matrix element of the NC between the initial and 
final hadronic states, and it includes all the nuclear effects. Introducing the variable y = q^ /E^ with E the incident 
neutrino energy, we can write 

da GP' t:^ f 1^ \ I 2 daL 2 da^ ^ das , dA 



dq'^dydtd(j}k^q IGtt^ \ \q\'^) ^ dtdipk^q^ dtd(j)k^q^ dtd4>k^q^ dtdipk^q^ ' ^ ^ 

where 

0^ 9' E + E'± \q\ 
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Note that 



das 1 dfjs dcjji 1 dan da^ 1 da 



L 



dtd(f>k^q 2TT dt ' dtd(j)k^q 27r ' dtdipk^q 27r 



(7) 



as neither Tioo: ''^zz nor the combinations TYqz + ''^zO: T^xx + ^yaj Wa;y ~ T^yx can depend on 0fe^g. The only quantity 
that depends on (fik^q is ^^^^ — . The latter is not a proper differential cross section as it can take on positive as well 

as negative values. Besides, it cancels upon integration on (t>k^q (see below). Note also that for ~ Q only the 0-5 
term contributes. 
The tensor H^"" , 

^f^"^ = j (8) 

only depends on g^, = (A^,0) (the four vector of the initial nuclear state, with A1 the nucleus mass) and t. Due 
to the tensorial character of H and the fact that, being the transferred momentum q aligned with the z— axis, it is 
invariant under rotations around the z— axis, one can prove 

j^xx _ jjyy jjxz _ jjzx _ jjyz _ jjzy _ jjOx _ jjOy _ jjxO _ jjyO _ q jjxy _ jjyx ^g^ 

Thus, integrating in (j)k^q one obtains for the ^^i^y^i differential cross section (32j 

TP f q^ \ r 2d(^L ^ 2d^R , o '^^s. i.^. 
;E K [u —— + V —— + 2uv——) , (10) 



where 



dq^dydt IGtt^ V I? IV dt dt dt 



^ ^ -4-(|gTi^oo + 9"l'fl(i^oz + i^zo) + (0'i^zz) , 

dt q^ K 

daR _ TT 
dt K 



{H,,-iH.,y). (11) 



As stated above, at = only as contributes, and given that in this case q'^ — \q\, one finds that the cross section 
goes as 

{\qfHoo + I (Hoz + Hm) + {qy H,,) = q^q'^H^, . (12) 

In other words, in the q^ = Q limit, the lepton tensor of Eq. ([3]), turns out to be proportional to q^q"^ and thus, one is 
left to compute the matrix element of the divergence of the hadronic current. Since the vector NC is conserved, only 
the divergence of the axial part contributes to the cross section. PCAC can then be used to express the divergence of 
the axial current in terms of the pion field operator ^{x), 

d^A%cix)^2Ummx), (13) 
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with /tt ~ 92.4 MeV, the pion decay constant. Treating the nucleus as an elementary particle, one can write 



{N-gs7T°ik^)\q^A';fc\^gs).o = ~2iUT{jVgs7r°ik^) ^ ^°(g)M 



(14) 



where T(/ <— i) is the transition amplitude between the initial hadron state plus a pion of four momentum q^, and 
the final hadronic state. Using this relation, one obtains 



das 
dt 



= - 4- 



. E-ir 



dt 



(15) 



and then, neglecting the nucleus recoil (q'^ = E.^), one can further write (Adlcr's PC AC formula (ssj ) 



da 



a^fl Euv da{Tr°Afg, ^ ttX.) 



dq^dydt 

In the 5^ = limit that we are using 
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Euv 



(16) 



(17) 



and thus 



da 



dx dy dt 



G^ME 



dt 



g2=0,E^=gO 



(18) 



with X = —q^jlMtf' and M the nucleon mass. This latter form, Eq. HH]), was adopted in the original RS model p^ . 
Recently, Berger and Sehgal [ll] have proposed to use Eq. ^6]) instead, for finite . To our understanding, both 
choices are somehow arbitrary. As shown in the derivation, there are several factors of the type which do not 

affect the = calculation, and that might lead to different corrections far from that limit. 

An important remark is in order here. Adler's PCAC formula relates the neutrino induced cross section to the 
off-shell (0 = 7^ = fc^) elastic pion-nucleus one (see Eqs. (fT5)) . (fT6|) . and p8)) ). It is tempting to neglect off-shell 
effects and approximate the latter cross section by the experimental one, or by any realistic model for it. This can 
be only strictly correct in the case of a pointlike nucleus. This was firstly pointed out by J.S. Bell shortly after 
Adler proposed his PCAC formula. Because of absorption and inelastic collisions, physical pions do not penetrate 
into the interior of heavy nuclei, and thus the pion-nucleus elastic cross section is only sensitive to the nuclear surface 
(roughly it scales as ^ with A the nucleus mass number). But neutrinos penetrate to all parts of nuclei, being then 
sensitive to the whole nuclear volume; for them cross sections scale as A [3J|. Thus Bell expected, that although the 
on-shell (t^o_/v^ cross section is surface-like, the off-shell (g^ = 0) one turned out to be volume-like. That clearly hints 
at a non-trivial off-shell behaviour for a^^aj^. Similar arguments were raised and discussed by Bernabeu, Ericson and 
Jarlskog in the context of muon capture in nuclei [35] . 

With all these caveats, there are cases where approximating the off-shell pion-nucleus cross section by the on-shell 
one in Eqs. (fTS]) . p6|) . and p8|) could be accurate and easily obtained in any good model for pion nuclear scattering. 
For instance, i) long wavelength limit: the pion wavelength is larger than the nucleus and it probes the whole nuclear 
volume (this in the case of the nuclear beta decay), or ii) short wavelength limit: here the pion interaction is weak 
and there is little multiple scattering. Practically we have a collection of independent scatterers (nucleons). 

This problem has not been taken into account in the recent works of Refs. [Tgj . where it is proposed to use 
the experimental elastic pion-nucleus cross section in the Adler's PCAC formula at energies where the process is 
dominated by the weak excitation of the A(1232) resonance and its subsequent decay into a ttTV pair. For resonant 
energies, the pion wavelength is such that it renders doubtful the assumption of a pointlike nucleus, while the pion- 
nucleon interaction is sufhciently strong to expect that surface effects due to the distortion of the incoming pion waves, 
included in the experimental elastic pion-nucleus cross section, might induce inaccuracies in the computation of the 
weak process. 

In order to remove such an effect from the physical a^aj^ one should rely on a distortion model, for instance the 
one described below in Subsect. IIVBI The original model of Ref. [Ij] deals with this volume- surface problem when 
< a: >, the average path length of the pions, is calculated by choosing the path length traversed by pions uniformly 
produced in the nucleus instead of the path length of pions scattered on the nucleus. 
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III. REVIEW OF THE REIN-SEHGAL MODEL 

As discussed above, the Rein-Sehgal model [T^j uses the Adler's PCAC formula [s^ and approximates the coherent 
7r° production differential cross section in the laboratory frame by means of Eq. ([T5|) . Besides the neglect of the 
off-shell (g^ 7^ mj) effects in the right hand side of Eq. (fT^ . this expression already involves a few approximations, 
that will not be discussed in this work. For instance, the final nucleus recoil energy is also neglected, which allows to 
approximate q° by the pion energy E.^. 

In Ref . , the expression of Eq. ()18p for neutrino coherent pion production is continued to non-forward lepton 
directions {q^ ^ 0) by attaching a propagator term 1/(1 — /m\)'^ with ttla ~ 1 GeV. This does not mean that the 
dependence is merely given by the added propagator term, and in fact it is much more pronounced than that induced 
by this factor. This is because large (and negative) q^ values are suppressed by the elastic pion-nucleus differential 
cross section that strongly favours t — Q. Within the = E'jr approximation, a zero value for t would imply = m^, 
a kinematical point that cannot be reached since q^ is space-like. Thus, the lowest values accessible for t must come 
from the w region. 

Once q^ ^ there are other finite contributions stemming from , cfr and additional pieces of as fully disregarded 
in the RS model. These contributions come from both the axial and the vector part of the NC. These further 
approximations have been, somehow, justified because of the strongly forward peaked character of the process that 
implies that only quite small values of q^ will be relevant. This assumption will be discussed later in more detail. In 
particular, as the contribution of the vector current is neglected, the model predicts equal neutrino and antineutrino 
cross sections. 

Next, in Ref. fl5|, the pion-nucleus cross section, with the caveats mentioned above, is expressed in terms of the 
pion-nucleon one and then one obtains 

da G^ME , 1 /,^,m2^ daiT:^N ^-k^N) \ 



dxdydt 7r2 •"^^ (l-q2/l GeV2)2 V dt 

where F_/\{t) is the nuclear form factor which can be calculated as F_^{t) — J d^r e^^^''^"^'^ {Ppif) + Pn(^)}, with the 
density Pp(„) normalized to the number of protons (neutrons). Finally, according to [l^ , _Fabs is a independent 
attenuation factor that takes into account effects of the outgoing pion absorption in the nucleus^. 

At this point, some new approximations have been implemented to estimate the pion-nucleus cross section. First, 
the pion-nucleon cross section is evaluated at t = 0, namely, in the forward direction. This can be justified if the 
nuclear form factor is sufficiently forward peaked. The larger the pion energy and the heavier the nucleus, the better 
this approximation becomes. In the original paper p^ . Eq. (I19p was employed for a medium size nucleus, aluminum, 
and neutrino energies above 2 GeV, for which the relevant pion energies are quite high. However, as we will show 
below and it was already pointed out in (25j . for neutrino energies below 1 GeV and lighter nuclei, like carbon or 
oxygen, the nuclear form factor is not enough forward peaked to render the finite t— dependence of the pion-nucleon 
cross section negligible, and even in the forward direction the t value is not close enough to zero. Second, the distortion 
factor -Fabs is an oversimplification since in any realistic scattering model this factor should depend on t. 



IV. IMPROVEMENTS ON THE RS MODEL 



There exist sophisticated microscopical calculations of the neutrino coherent vr" production off nuclei based on the 
dominance of the A(1232) resonance [l^, [2l|, [l^, [H, [l^, H^- However, these are difficult to extend for the higher 
neutrino energies present in the current experiments and to implement in Monte Carlo algorithms. Actually, one of 
the virtues of the RS model is its simplicity, that allows its use for different nuclei and pion energies. This is one of 
the reasons why the RS model has been widely used. 

We propose here a minimal extension of the model in which we implement two main corrections. Firstly, we will 
take into account the dependence of the pion-nucleon cross section when computing the pion-nucleus scattering, 
and secondly we will work out a more realistic description of the outgoing pion distortion^. 



^ As defined in T^, -Fabs only removes from the flux pions that undergo inelastic collisions but no true pion absorption is actually 

accounted for by means of this factor [25ll . 
^ The pion-nucleus cross section can not be factorized, in addition to the vr N elementary cross section, as the product of a nuclear form 

factor (defined as the Fourier transform of the density) and a distortion factor. This is in contrast with purely clectroweak processes 

not involving strong final state interactions. 
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A. f-dependence of the pion-nucleon cross section 



The t — approximation in the right hand side of Eq. (jl9p is unnecessary and can be easily removed, and one 
obtains 



da 



dx dy dt 



dt 



(20) 



where cr„s/ is the non spin-flip part of the pion nucleon cross section because spin-flip processes do not contribute to 
the elastic pion nucleus cross section. Given that for t = the spin-flip cross section is identically zero this constraint 
was unnecessary in Ref. [l^l- The non spin-flip cross section can be evaluated in terms of the phase shifts ((5/^) and 
inelasticities i), that we take from Ref. [sHi (J, ( and / stand respectively for total angular momentum, orbital 
angular momentum and isospin of the pion-nucleon pair). For an incoming neutral pion of laboratory energy E-m the 
non spin-flip cross section reads 



dcr- 



nsf 



with s 



M2 



dt 
a^{s,t) 



2ai 



^l^i'lc.m. 



£-1/2 



Pt (cos ( 



[M + r, 



cosf 
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cm. 



(21) 



(M- 



/4s, the pion momentum squared and P( 



the Legendre polynomials. For the sake of completeness we also give here the expression for the inelastic pion-nucleon 
cross section, Omei 



lis) 



I I c.r 



lm(|ai(s,t = 0) + ia^(s,t = 0)) -47r^ [(^+l)|a 



,^+l/2|2 



e\a 



'-1/2|2 



(22) 



with a/(s) = ^2 f^^^ ei^) + fi/2 ei^)) '^'^ discuss the importance of the t dependence of the cross section we plot 
in Fig. [T] the ratio W{t) /W{0) for carbon and iron at two different pion energies . The function W{t) is defined as 



W{t) = \FAt)\ 



2dansf{T^"N 



dt 



(23) 



and for the RS model 



rfgns/(7r°jV->7r°jV) 



dt 



is taken at t = 0. To allow for a better comparison with the RS model we 



have used the same nuclear form factor as in Ref. given by Fj^{t) = Ae''*^^, with b = R^A^^^ and A the number 
of nucleons. As in Ref. [l^ we take Ro = 1 fm. As can be appreciated in the figure, the extra dependence of the 
non-spin flip pion-nucleon cross section on t, neglected in Ref. |12l |. produces drastic changes in a light nucleus like ^^C 
for low energy pions. We have chosen iJ^ = 250 MeV because it corresponds approximately to the maximum of the 
7r° yield in the MiniBooNE experiment [1^ . The corrected distribution is much narrower and leads to a signiflcantly 
smaller area. Effects become less important for increasing pion energies and for heavier nuclei thanks to the more 
pronounced behavior of the nuclear form-factor. In the original work of Ref. [12], aluminum and neutrino energies, at 
least of 2 GeV, were considered. For 2 GeV neutrinos, the pion spectrum will peak around 1 GeV, and for such pion 
energies the corrections to the RS formula turn out to be moderately small. But as seen, they can be relevant for the 
low energy pions produced in low energy neutrino processes. 



B. Pion distortion 



Next, we pay attention to the distortion of the outgoing pion. We use an eikonal approximation, in the spirit of 
the original RS model, and taking q in the z direction, we replace 



|i^^(i)l'^^abs 



dVe'('?-'^")->(r) 



P 



abs 



(24) 
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FIG. 1: W{t)/W{0) ratio for carbon and iron at two pion energies. 
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with p{r) the nuclear density, b = \/x^ + y^ the impact parameter, = \q \ — q - k-^/lql, pt 
and Jo the Bessel function. The eikonal distortion factor T{b, z) is defined as 



z) = exp {--(Tinei J dz' p y^Jb"^ + z'2j 



(25) 
(26) 



{q- KY/q- 



(27) 



If one takes r(6, z) independent of b and z, F^^^°^ {\q\, |, tf- fc^r) depends only on pI~^Pt — ~^ ^^"^ reduces to the 
structure originally proposed in Ref. [l^, where all t dependence is ascribed to the nuclear form factor Fj\^{t). Indeed, 
in the original RS model some volumetric average for T(b, z) is assumed. Neglecting the i— dependence, inherited from 
the impact parameter dependence of the distortion of the pion waves leads to flatter angular distributions as we will 
show. 



V. PION NUCLEUS SCATTERING 



In this section we will confront with data the predictions of the two previous approaches when extended to elastic 
pion nucleus scattering. By studying low energy pion-nucleus scattering we magnify the effect of the distortion factor 
and in this way we expect to gain some insight into its possible deficiencies when applied to low energy neutrino- 
induced coherent pion production. 

Following Ref. [13], one would write for the pion nucleus clastic cross section 



dt 



da{n°N 77° N) 



dt 



(28) 
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where we have replaced the original absorption factor i^abs, appropriate for neutrino induced pion production, by 
that also takes into account the initial pion distortion. We have done this following Ref. [l^l- Namely, the original 
absorption factor Fabs given in Ref. reads: 



(29) 



where < x > is the average path length traversed by a tt'^ produced homogeneously in the nuclear volume by the 
neutrino, and A^^ = (JmeiP with (Tinei the inelastic Tr^-nucleon cross section. The average < x > is found to be 3i?/4 
assuming a hard sphere density of radius R for the nucleus. Thus, taking R — Ri^A^ , one finds 



abs 



T ( 



(30) 



To properly compare with pion nucleus data one should also consider the distortion of the incoming pion. In this 
case, after averaging over the impact parameter one obtains < x >= 4i?/3 and thus we replace i^abs by F^-^^ where 



ll/3 



abs 



(31) 



Using the original Fabs instead of F^^^ amounts to a change of scale but does not modify the angular shape of the cross 
section. Finally, we obtain atnei and '^-(-°^--°^) fr om the SAID partial wave analysis [36|; detailed expressions 

were given before. Before discussing the results, we would like to point out that the calculation is for 7r"-nucleus 
scattering while data correspond to charged pions. However, we have only selected data from isoscalar nuclei and 
except for Coulomb effects (relevant at small angles) the cross section should be the same. 

We would like to stress the fact that the authors of Ref. [l2| did not actually make a model for low and medium 
energy pion-nucleus scattering. The extension to those energies that we are showing here should certainly not be 
attributed to the RS model. 

What we see in Fig. [5] (solid lines), for all nuclei and pion energies shown, is that both the size and the angular 
dependence predicted by this model derived from the RS approach strongly disagree with data. We have explored the 
available data for other nuclei and pion energies and the outcome is similar. Within this model, the size of the elastic 
cross section could be largely affected by even small variations of i?o (see Eq. ()31|)). and moreover it is not very clear 
why one should use cFinei in the absorption factor, as we discuss below. 

In the case of coherent pion production induced by neutrinos one expects distortion effects, accounted for by Fabs, 
to be less relevant and besides one should bear in mind that experimental analyses do often adjust the size of the cross 
section with a free parameter. However, those neutrino analyses always rely on a theoretical model for the outgoing 
pion angular distribution. As shown in the figure, the RS derived model leads for low pion energies to 7r-nucleus 
differential cross sections much flatter than experiment for all nuclei and one would expect this will also be the case for 
coherent pion production induced by low energy neutrinos. From this latter perspective, we tentatively conclude that 
the RS model for neutrino-induced coherent pion production, widely used in the literature, might induce important 
uncertainties in the analysis of the neutrino oscillation experiments for low energy incident neutrinos. 

We discuss now the results that one obtains with the modifications to the RS model introduced in section [TVl The 
pion-nucleus scattering would now read 



where 



dt 



dt 



dVe'('--^")-V(f)r'^(6) 



r^b) = exp(^-ia,„e( dz'p ( + z'^) 



(32) 

(33) 
(34) 



The profile F'^ does not depend on z since both the incoming and outgoing pions should now be distorted. 

We show our predictions (Eq. (|32p ) in Fig.[2]with dashed- lines. As can be seen in the plots the two simple corrections 
(i— dependence of the pion-nucleon cross section and a more realistic description of the pion distortion) lead to a better 
description of the angular dependence, improving enormously previous RS based model results. 

Finally, we would like to point out that it is not clear why one should use atnei in the distortion profile. Indeed, 
within the eikonal approximation one should use the imaginary part of a realistic pion-nucleus optical potential [23 | . 
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FIG. 2: Pion nucleus elastic cross section for several nuclei and energies. Solid line: RS based model. Dashed line: Improved 
model of Eq. (I32p . For comparison in the right top panel, we also display results (dashed dotted green line) obtained from the 
solution of the Klein-Gordon equation with a microscopic potential calculated within the A— hole model jSTJ], as employed in 
Refs. [231 . 113, H^. To compute our results, we have used modified harmonic oscillator [p(r-) = po{l + a{r/R)'^) exp( — (r/f?)^)] 
densities for carbon and oxygen and a two parameter Fermi distribution [p{r) = Po/(l + exp((r — c)/ao))] for calcium. Density 
parameters: R = 1.692 fm, a = 1.082 and R = 1.833 fm, a = 1.544 for ^^C and ^'^O; c = 3.51 fm and ao = 0.563 fm for *°Ca. 
Data are from Refs. '3S] C^C), ^ (^''O), H (*°Ca at T„ = 230 MeV) and 5l| ("'"Ca at T„ = 673 MeV). 



This would eliminate those pions that are absorbed, undergo quasielastic processes or suffer inelastic reactions, tak- 
ing into account Pauli blocking, Fermi motion and other many-body effects to evaluate the corresponding reaction 
probabilities. Such kind of sophisticated optical potential is only available for low energy pions, up to the region of 
the A(1232) resonance (see for instance Refs. [33, H^, H^), and it has been used in the microscopical approaches of 
Refs. i2l[2a,[24, 2 5]. Its extension to higher energies is a highly non-trivial task. Here, we are not so much concerned 
with the size of the cross sections but rather with the outgoing pion angular and energy dependence. For this latter 
purpose, it might be sufficient to use ainei to compute the distortion profile. It is true that in this way, we do neither 
account for pion absorption nor for pion quasielastic distortion (induced by pion-nucleon elastic scattering). However, 
this is partially compensated since the use of cTmei leads to a larger distortion than it would be expected, if the inelastic 
processes were considered in the nuclear medium and Fermi statistics was taken into account. The same discussion 
applies to the profile function, T{b,z), relevant for the neutrino coherent pion production reaction, and defined in 
Eq. 



VI. NEUTRINO-INDUCED COHERENT tt^ PRODUCTION RESULTS 



Once our framework has been satisfactorily tested against elastic pion-nucleus differential cross section data, we 
give here results for neutrino-induced coherent tt^ production. We work in the neutrino-nucleus LAB frame. For the 
X, y, t differential coherent tt" production cross section we use 



da 



dx dy dt 



G'^ME , i/[l-|cosl 



(l-g2/lGeV^) 



2^2' 



\K 1,9 - K)\ 



^dansfin°N^TT''N) 



dt 



E„=q° 



(35) 



with H[...] the step function and 9 the center of mass pion angle in the free pion-nucleon elastic reaction (see Eq. (PT|) ). 
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FIG. 3: Plane Wave (no pion distortion) distributions obtained from models I and II, together with those deduced from the 
RS and RS* models. In addition for E=0.5 GeV in carbon, we also display results from the microscopic approach described in 
the text. Filled (open) squares stand for the calculation with the full V-A (only with the contribution) A^'A current. 



That defines our model I. We will also consider a second model (II), where some kinematical corrections, recently 
proposed by Berger and Sehgal [1^, are incorporated. Those corrections, and their degree of arbitrariness, were 
already discussed in Section |TT] and they amount to replace in Eq. (|35p 

We also incorporate these corrections in the original RS model performing the above replacement in Eq. (|19p . In what 
follows, we will refer to this latter model as RS*. Although we will look at the effect of the kinematical corrections 
proposed in Ref. [ll], we will not present results obtained from the full approach followed in that reference. The 
reason being that in Ref. [l^ the use of the experimental pion-nucleus elastic cross section is advocated. This is not 
fully correct because of the strong distortion of the incoming pion, implicit in the experimental cross section data, 
that should not be taken into account in neutrino-induced coherent pion production. 



A. g^— distributions 

We readily obtain da/dq^ from Eq. (j35p by integrating over y and t and performing the change of variables x ^ . 
In Fig. [3l we show results obtained without distortion of the final pion. Eliminating distortion allows us to check, by 
comparison with a microscopical calculation, the validity of the t — approximation used in the RS and RS* models. 
With that aim, in carbon and for E=0.5 GeV, we also display results from the microscopic model of Ref. [1^. At 
this energy, neutrino-induced coherent pion production is dominated by the excitation of the A(1232) resonance and 
its subsequent decay into a tt'^A^ pair. Thus, for simplicity, we have neglected the tiny effects of the non-resonant 
background and the crossed— A term included in Ref. [2^. Besides, as all discussed models are based upon the 
free pion nucleon cross section, nuclear medium effects that modify the A propagator have not been included. We 
consider two cases, one that takes into account the full A^A weak current (V-A Ap), and another that only includes 
the dominant axial term proportional to the form factor^ (Only Ap). PCAC based models rely on Eq. 



* See for instance Eq. (1) of Ref. | 44|| for a form factor decomposition of the A^A weak current. 
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which relates, at = 0, a weak matrix element with a purely hadronic one. In the case of the NA transition, the only 
term of the divergence of the axial current that survives at — is that proportional to the form factor C^. Thus, 
to make meaningful the comparison between the PCAC based models examined here and the microscopic approach 
including explicitly a A, requires fixing C^(0) to the value predicted by the non-diagonal Goldberger-Treiman relation 

(C^(0) = y|/7r£; = 1-2, with the nNA coupling /* = 2.2 fixed to the A width). We see that I and II calculations 
are in reasonable agreement with the microscopical model, whereas the RS and RS* ones predict much larger and 
wider differential cross sections. As mentioned above, this is mainly due to the t — approximation assumed in those 
models^. The effects due to the t = approximation decrease with the neutrino energy and atomic number, but they 
are still important for E=2.1 GeV in a medium sized nucleus like calcium. At q^ — 0, and in the absence of distortion, 
our model should essentially match the microscopical calculation (as it was discussed in Sect. IVB-1 of Rcf. 25]), 
and indeed the two models nicely agree at this kinematical point. 

Far from q^ = 0, PCAC based models just consider the axial part of the current, that given the tiny effect of the 
vector part in the microscopical model (as can be seen by comparing "V-A A^" and "Only Ap" microscopic 
results), it turns out to be an excellent approximation here. Our model II describes the "Only Ap" results better 
than model I, which make us conclude that the kinematical corrections recently proposed by Berger and Sehgal [l3| 
are sound. 

The small differences at finite q^ can be due to some discrepancies between the experimental pion-nucleon cross 
section and that deduced from the A mechanism alone (with /* fixed from the A width), to offshellness {q^ m^) 
effects, or to different treatments of the nucleon dynamics inside the nuclear medium. Besides, there are other effects 
that vanish at = and that are present in this case. For instance, in the microscopical model the as structure 
function (Eq. pT|) ) is not proportional to q^^q'^H^i, and hence as cannot be exactly related to the pion- nucleus cross 
section for ^ values. Similarly, corrections due to the an, a^ structure functions in Eq. IjlOp are not taken 
into account within the PCAC approximation. Other sources of discrepancies come from the adopted q^ behaviour 
of C^{q^) [1^, which gives rise to a faster g^— decrease than that provided by the propagator term 1/(1 — g^/m^)^ 
included in the PCAC based models. Nonetheless, at least for this observable, the effect of all these differences is 
minor and model II gives a good approximation to the microscopical calculation^. 

Now we look at the results with distortion of the final pion that we show in Fig. [D In this case, for comparison, we 
use the model based on the microscopic approach developed in Refs. [13, including both distortion and the nuclear 
medium effects affecting the A(1232). Again, and for simplicity, we have neglected the tiny effect of the non-resonant 
background and the crossed— A term included in Ref. j25j|. 



^ The used nuclear form-factors are also different (see text after Eqs. II19I I and l(23j), but their effect can not account for the large differences 
observed. As a matter of example in ^^C and for E=0.5 GeV, the plane wave RS model gives an integrated cross section (in units of 
10~^* cm^) of 0.206, while our model I predicts a value of 0.070 in the same units. If we use model I with the nuclear form factor 
proposed in the original RS model, we find an integrated cross section of 0.089 10~^^ cm'^. 

® A word of caution must be said here. The basic PCAC formula in Eq. I I14I I might suffer from corrections. For instance, in the A region, 
a value for C^{0) significantly smaller (0.87) than 1.2 was used in Ref. 25]. This reduced coupling would lead to cross sections around 
a factor of two [(1.2/0.87)^] smaller than those presented in Fig. |3] This violation of the off-diagonal Goldberger-Treiman relation by 
about 30% was proposed in Ref. 9] after fitting the Argonne bubble chamber -I- p — » ^~7r+p cross section data [Toll . including a 
non-resonant background. Non-resonant background terms, though important at the nucleon level, turn out to be negligible in the 
neutrino coherent pion process |25[ . The lattice QCD results shown in Figure 4 of Ref. [45ll might also support a value for C^(0) smaller 

than \/^/vr^— ■ This is indeed a serious problem that deserves a joined theoretical and experimental effort to sort it out. 
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FIG. 4; Full distributions (including distortion of the outgoing pion). For E=0.5 GeV in carbon we also display results from 
the microscopic approach of Ref . [H with C^(0) = 1.2. Filled (open) squares stand for the calculation with the full V-A (only 
with the contribution) A''A current. 



Distortion reduction is significantly larger in the RS and RS* models than in our scheme. That diminishes the 
differences in the full calculation in carbon, and explains the reverse pattern observed in calcium. Note, however, that 
the pion distortion in the RS and RS* schemes is too strong. Indeed, it leads to forward pion-nucleus elastic cross 
sections smaller than both data and those deduced within our framework, as can be appreciated in Fig. [21 In spite of 
that, both the RS and RS* models still overestimate the microscopical calculation in carbon by a large factor. 

On the other hand, the situation is somehow puzzling when we compare microscopical results with those of model 
I or II. As discussed above, within the PW approximation both types of approach nicely agree, but distorted results 
notably disagree. This is despite the fact that the Klein-Gordon A— hole model, in which the microscopical results 
depicted in Fig. [5] are based, leads to elastic pion nucleus cross sections of similar quality as those deduced from the 
eikonal framework developed here, and that our model I and II use (see for instance the comparison performed in the 
right top panel of Fig. [2] for a typical pion energy). This contradiction deserves some discussion, and we believe it 
is related to the inadequacy of the use of eikonal distortion for low energy pions and to the modifications of the A 
resonance properties in the nuclear medium. 

Describing correctly the elastic pion-nucleus cross section is a necessary condition that should satisfy any model 
aiming at a proper description of the coherent pion production process induced by neutrinos. However, that constraint 
is not sufficient, and the situation in Fig. |3]is an example of that. 

One of the major sources of differences between our model of Eq. ([5^ and the microscopic calculation is due to 
in- medium modifications of the free pion-nucleon cross section in the latter case (see for instance Ref. [1^). Within 
the eikonal model assumed here to construct the elastic pion-nucleus cross section, there are two competing factors: 
the pion-nucleon elastic non spin flip (Jnsf cross section and the pion distortion, controlled by <J\nc\ ■ As a consequence, 
the pion-nucleus scattering amplitude does not linearly depend on the corresponding pion-nucleon one. That explains 
why different approaches for the in medium pion-nucleon amplitude could lead to similar pion-nucleus cross sections. 
For instance, let us suppose that in medium modifications reduce (enhance) the pion-nucleon cross sections. We 
have that pion distortion effects become less (more) important with decreasing (increasing) CTinoi and this reduced 
(enhanced) suppression can compensate a decrease (an enhancement) in Unsf- The outcome is that you can predict 
similar pion-nucleus cross sections starting from quite different pion-nucleon cross sections. It seems that induced 
differences in the pion distortion compensate the changes due to the medium in the elastic pion-nucleon cross section 
in pion-nucleus scattering, but not in neutrino induced coherent pion production where distortion becomes smaller 
and affects only the outgoing pion. 

For low energy pions, of interest in MiniBooNE and T2K, we believe that the distorted results from the microscopical 
models of Refs. [2J, [2^ are more realistic that those based on the eikonal approximation presented here, where the 
pion-nucleon interaction is not modified in the medium and a simple model based on Uinei is used to distort the 
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outgoing pion waves. The framework of Ref. [24], despite the use of the eikonal approximation, greatly overcomes 
these shortcomings, since there, the modification of the A properties in the medium are taken into account, and the 
imaginary part of a reahstic pion-nucleus optical potential is used to distort the outgoing pion. The model becomes 
as complicated as that of Refs. [13, HE], as difficult as this latter one to extrapolate at higher pion energies, and it is 
still less reliable, since multiple scatterin g is not taken into account within the eikonal approximation as accurately 
as by solving the Klein Gordon equation [24 . [25l | . 



B. Pion distributions 



For NC driven processes, the distribution can not be easily measured because of the obvious difficulty in detecting 
the outgoing neutrino or the nucleus recoil. Thus, it is of great interest to examine pion distributions. 



1. Pion energy distributions 



In Fig. O we show the outgoing pion spectrum from oxygen, with and without distortion effects, predicted by the 
different models considered in the previous subsection. The incoming neutrino energy is 0.65 GeV. The integrated 
cross sections are compiled in Table HI 
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FIG. 5: Distorted (left panel) and plane wave (right panel) da/dET^ results in oxygen obtained from our models I and II, 
together with those deduced from the original RS and the RS* models. The incoming neutrino energy is 0.65 GeV. In addition, 
we also display results from the microscopic approach of Ref. [25| with C^(0) = 1.2. Filled (open) squares stand for the 
calculation with the full V-A (only with the contribution) A'^A current. 

We see once more that, for the case without distortion, our model II has a reasonable agreement with the "V-A 
Ai?" and "Only C^Ap" microscopic calculations. Not only the shapes, which peak at around E-^ = 320 MeV, but 
also the integrated cross sections are similar (see the table). Again the RS and RS* models give much larger results 
and the position of the peak is slightly displaced towards lower energies. 

In the distortion case we find totally different behaviours. Results obtained with our models I and II show a shape 
similar to the case without distortion, while the RS and RS* calculations show a double peak structure. In the 
microscopical calculation a single peak appears, but now at a lower pion energy i?^ « 260 MeV. This change in the 
position of the peak in the microscopical results can be traced to the modification of the properties of the A in the 
medium and to nonlocal effects in the amplitude (see for instance Fig. 2a in Ref. [l^). These non- localities reflect the 
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Microscopical 








This work I 


This work II 


V-A A Only A 


RS 


RS* 


PW 


0.152 


0.127 


0.117 0.106 


0.389 


0.305 


DW 


0.072 


0.060 


0.037 0.033 


0.100 


0.078 



TABLE I: Integrated cross sections (in units of 10 cm^) corresponding to the different curves displayed in Fig[5] PW and 
DW stand for plane wave and distorted results, respectively. 



fact that the pion three- momentum is only weU defined asymptotically when the pion- nucleus potential vanishes [25l | . 
These efi^ects are neither present in any PCAC based model where one uses the free-space pion nucleon cross sections, 
nor in any microscopic calculation with plane waves. 

At low and intermediate neutrino energies, once the distortion is considered, we must again conclude that the 
microscopic model predictions are most reliable than those deduced from the PCAC based model examined here. The 
original RS predictions have little reliability and the hope is that at higher energies, where we expect nuclear medium 
modifications of the elementary pion-nucleon interaction to be less relevant, the model I and II could become better 
suited. 



2. Pion angular distributions 



The starting point of the PCAC based models is the x, y, t differential cross section. Firstly, we would like to stress 
that for 7^ 0, the knowledge of ^^^"^^^ is not sufficient to compute the angular distribution of the outgoing pion with 
respect to the direction of the incoming neutrino. Let us take here as z-axis the direction of the incoming neutrino. 
The pion and transferred momenta arc given by 

kjr = |fc7r| (sinfl^ COS 0^, sin^^ sin0^, COS 07r) , (37) 
q= f-|fc'|sin6'',0,|fc| - |fc'|cos6l') , (38) 



with 9' the outgoing neutrino scattering angle, the angle between the incoming neutrino and the outgoing pion 
and the azimuthal pion angle in this frame. By means of a rotation that takes q along the positive 2— axis one can 
easily obtain 

- - I A:' I I fc I — I fc' I cos 9' 
cos 9t^ — k ■ — sin 9' sin 9^ „ cos „ H — cos 9^ „ , (39) 

where 9k^q and (t>k^q are the pion polar and azimuthal angles in that frame where the positive z— axis is taken in the 
direction of q. 

The incoming neutrino energy and the variables x and y determine |fc'|, \q\ and 9', while, within the = q^ 

approximation, t fixes 9k^q [t = —q^ — fc^ + 2\q Wk^rl cos9k^q ], remaining (j)k„q undetermined. Thus, as stated before, 
the knowledge of da/dx dy dt is not sufficient to compute the angular distribution of the outgoing pion with respect to 
the direction of the incoming neutrino, and it would be necessary to know dtd<j>k — ' ^'^^^ ^'^^ 1^ ~ ^' ^k^i ~ ^'^^ 
since 9' = (both the outgoing neutrino and the momentum transfer go along the incoming neutrino direction) and 
, , J , r. determines , , '^'^ — 3- 1 , „. 

ax ay at ] —0 ax ay a cos \q^—0 

In PCAC based models it is assumed, 

da 1 da 1 f da 

........ ' (40) 



which leads to [12 



dx dy dt d(f>k^q 2tt dx dy dt 2n J " dxdydtaqjk^q 



da E-E^\k^\\q\ (AA\ 

-ii [1 - I cost/^IJ — — — - , (41) 



dEj^d-q d cos 9' dcfi-K ME^ tt dx dy dt 

where = i?^(l — cos6'7r) is the variable proposed by the MiniBooNE Collaboration in its recent analysis of coherent 
7r° production of Ref. ,2^] . 
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FIG. 6: Undistorted da/drjdq^ (top panels) and da/drj (bottom panel) differential cross sections obtained from models I, II, 
RS and RS* and the microscopic approach of Ref. [25|] with C7^(0) = 1.2. Filled (open) squares stand for the calculation with 
the full V-A (only with the contribution) A*'A current. In the top-right and bottom panels we also show results (triangles) 
for the "V-A Af" microscopic calculation including only the as contribution. 



For non- vanishing values, Eq. is incorrect, and therefore Eq. (HT1) is wrong as well. The problem arises 
because dtd(t>k — depends in general on (t>k^q, through the ^^^^ — term in Eq. When is zero that term 
does not contribute and thus Eqs. (|40|) and ((4T|) are correct only in this limit. 

Let us start by looking at undistorted results in Fig. [6l In the two upper panels we show ^j;^^ differential cross 
sections for two different q^ values, while in the lower panels we show q^ integrated distributions. The nucleus is 
carbon and the incident neutrino energy is i? = 1 GeV. For q^ = (top-left panel) we find a very nice agreement 
between our models I and II calculation and the microscopic model results for j^^- On the other hand, RS and RS* 

models give rise to significantly flatter distributions. In the top-right panel, we show results for q^ — —0.25 GeV. 
We find now a disagreement of the models I and II with the microscopic calculation. This disagreement goes over 
to the q^ integrated ^ differential cross section shown in the bottom panel. The microscopic calculation is much 
more peaked at = than either models I and II or the RS and RS* models, the latter two showing the flattest 
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distributions. The reason for the discrepancy between the microscopic model and the PCAC based ones is due to two 
approximations made in the latter models. The first approximation made in PCAC based models is the one encoded 
in Eq. (j40| which leads to Eq. ([4T|) instead of the correct expression 

,^ , , ,, = 2^^— ^Ifc^l IqliJfl - I cos6l^n -—r4^-r- (42) 

dE^drjd cos e'd<j)„ ME^ ' ' ^ ' ""^^ dx dy dtd(j)k^q ^ ' 

As mentioned above, the approximation in Eq. (f40|) amounts to neglect the ^^^^ — term in Eq. fl]), term that the 
microscopic calculation properly takes into account. 

Besides, PCAC based models only consider the as contribution in Eq. ^ while the microscopic one also takes into 
account the contributions from an and aL- To see the effects associated to these neglected terms we also show in the 
top-right and bottom panels the "V-A Aj?" microscopic calculation considering only the as contribution. What we 
see is that the differential cross sections become much flatter than before, being now in a nice agreement with our 
model II calculation. The fact that the RS and RS* models give the flattest distributions is explained by the further 
t = Q approximation used there. 

The conclusion is that in PCAC models a^^ aL and A structure functions in Eq. (U cannot be taken into account. 
The omission of these contributions leads to significantly flatter da/drj differential cross sections than those predicted 
by a microscopic calculation. This affects both our models I and II and the RS and RS* models. In the latter two 
cases the t = Q approximation enhances this behaviour. 

This is also true for distorted results that are shown in Fig. [71 Again the microscopic calculation is more peaked 
close to = 0, and the RS and RS* results show the flattest distributions due to the t — Q approximation. In any 
case, model RS* and even more model II represent a major improvement as compared to the original RS model for 
this observable. 

These findings have immediate consequences for the results published by the MiniBooNE Collaboration in its 
recent analysis of coherent tt'^ production of Ref. [l^, which rely on the RS model. The distribution for E-n{l — cos 9t^) 
given there should be definitely much narrower and much more peaked around zero, and thus it might improve the 
description of the first bin value in Fig. 3b of this reference. This was already suggested in Ref. [25j . 



VII. SUMMARY AND CONCLUSIONS 



We have critically reviewed the commonly used Rein-Sehgal model for NC neutrino coherent pion production [l2j . 
We have unambiguously pointed out the main deficiencies of this model, which induce important uncertainties for 
pions of relatively low energy, as those of relevance in the MiniBooNE and T2K experiments. Among others, the 
more relevant ones are: 

1. The t dependence of the coherent production is fully ascribed to the nuclear form factor Fj({t), while further 
and significant t— dependences induced by the pion-nucleon interaction are ignored (see Fig. [T]). 

The recent works of Refs. [l^, [l^ try to overcome this problem by using experimental information on the 
t— dependence of the elastic pion-nucleus cross section. Apart from the obvious limitation coming from the lack 
of experimental data for many pion energies and nuclei, this might not be appropriate either for energies where 
the process is dominated by the weak excitation of the A(1232) resonance. As discussed in Sect. HIl there might 
be a non-trivial off-shell behaviour for tr^oj^/ [s^. In a microscopic approach one would argue that because of 
the strong distortion of the incoming pion in the on-shell elastic pion nucleus process at energies in the A(1232) 
resonance region, one cannot directly relate its amplitude to the pion production induced by a weak current. 
That takes us naturally to the second caveat. 

2. The treatment of the outgoing pion distortion within the original RS model is quite poor, and it turns out to 
be certainly insufficient for resonant pions. An unquestionable evidence for that can be seen in Fig. [21 where 
it is shown that both the size and the angular dependence predicted by a model derived from the RS approach 
strongly disagree with the elastic pion-nucleus differential cross section data. 

3. Far from the = kinematical point, any PCAC based model, and in particular the RS one, cannot be used to 
determine the angular distribution of the outgoing pion with respect to the direction of the incoming neutrino 
(see right top and bottom panels of Fig. [5]). Terms that vanish at = 0, and that are not considered in PCAC 
based models, provide much more forward peaked outgoing pion-incoming neutrino angular distributions. 

PCAC models can only determine the distribution on the angle formed by the pion and the lepton transferred 
momentum, q. Experimentally, one can have access to this latter differential cross section in the case of CC 
driven processes, but not when the reaction takes place through the weak neutrino NC. 
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FIG. 7: Same as Fig. [6] including distortion. 



We address the first of these deficiencies, and we succeed to improve the original RS model by incorporating the 
i— dependence of the pion-nucleon cross section. 

However, there is not an easy solution for the other two caveats at low and intermediate pion energies. 

We have tried to improve the treatment of the distortion of the outgoing pion, while keeping the model still 
reasonably simple. Although we have managed to describe the pion-nucleus elastic cross section, we still find significant 
discre pan cies in the case of neutrino induced processes when we compare with the accurate microscopical model of 
Refs. |24l. [25t (see Figs. |3] and [5]). This illustrates a further interesting point: describing correctly the elastic pion- 
nucleus cross section is a necessary condition that should satisfy any model aiming at a proper description of the 
coherent pion production process induced by neutrinos, but however it is not a sufficient one. Discrepancies are due 
to the modification of the elementary processes when they take place inside the nuclear medium, and to the highly 
non-linear character of the strong driven processes. 

Altogether, we have only been able to quantify the unavoidable systematic error associated to these two last caveats 
for low and intermediate energy pions by comparing with the model of Refs. (2^ . [25j . 

All its limitations notwithstanding, the improved models I and II give a much better description of low and 
intermediate energies coherent pion production than the widely used RS model. Of course, as the neutrino energy 
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increases, the = kinematics becomes much more dominant, and on the other hand we expect nuclear medium 
modifications of the elementary pion-nucleon interaction to be less relevant. Under these circumstances, our improved 
models become even more appropriate. Nevertheless, we should point out the existence of a real problem: for neutrino 
energies in the region 1 to 2 GeV, there does not exist a reliable model to describe the coherent pion production process. 
This is because, these neutrino energies are not large enough for PCAC models to properly work, and they are certainly 
in the limit of applicability of microscopical models which only include the dynamics of the A(1232)— resonance. 
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